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1. Introduction 



The problem of counting topologically distinct knots remains a challenging one: see 
[1] for a review and for a report on recent advances. It was recently noticed [2] that 
combinatorial methods developed in quantum field theory, namely Feynman diagrams 
applied to matrix integrals, may provide a new way to count knots. In particular, the 
counting of alternating tangles, which had been achieved in [3], was reproduced. This 
counting, however, is not capable of discriminating between objects with different numbers 
of connected components: in the usual terminology, it gives the number of links rather 
than that of knots. In the present article, we reexamine this question and show how the 
introduction of a number n of possible "colors" for the knotted loops would solve this 
question, and how this may be formulated in terms of a matrix integral. For n = 2, this 
integral is equivalent to one recently studied in detail and computed in the framework 
of the random lattice 6- vertex model [4-5]. We thus carry out the explicit counting of 
alternating 2-color tangles: their generating function is the solution of coupled equations 
involving elliptic functions. We are able to give the 13 first terms of its expansion and its 
asymptotic behavior. We conclude with more conjectural considerations on the number of 
(2-color alternating) links. 



2. A matrix model for colored links 

We want to consider a model which describes alternating links in which each of the 
intertwined loops can have n different colors. This can be achieved via a large N x N 
matrix integral [2]. As a first attempt, let us consider the following integral 

*•»{«,,) = ffl dM.e^ H M « + i SU- W„M t ) (2 i} 

J 0=1 

over N x N hermitean matrices, and the corresponding "free energy" 

Such an integral has a g-series expansion ( "perturbative expansion" ) which admits a graph- 
ical representation in terms of Feynman diagrams. In the large N limit, only planar dia- 
grams survive: see for example [6] for a review. The integral (2.1) has an 0(n)-invariance 
where the M a form the vector representation of 0(n). The Feynman diagram expansion 
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of F(n, g) generates planar diagrams with four-legged vertices and colored edges such that 
colors cross each other at each vertex (Fig. 1). 

To each planar diagram we can associate an alternating link diagram (see e.g. [7], 
page 21) by following colored loops as they cross other loops and choosing alternatingly 
under- and over-crossings. This can be carried out in a consistent way throughout the 
whole diagram. 




Fig. 1: A planar Feynman diagram of (2.2) and the corresponding alternat- 
ing link diagram. 

We have thus generated alternating link diagrams with n colors. Diagrams with 
different numbers of connected components can now be distinguished by the n dependence 
of their weight: indeed, to a diagram with k connected components is associated a factor 
n k . We can therefore write 

oo 

F(n,g) = J2 F k(g)n k (2.3) 

k=i 

where Fk(g) is the sum over alternating link diagrams with exactly k intertwined loops. 
Note that if one can define F(n, g) for non-integer n and in particular in a neighborhood 
of 0, one has access to the individual contributions Fk(g) since they form the small n 
expansion of F(n,g). As usual in 0(n) vector models, one can perform (at least formally) 
an analytic continuation in n by a Hubbard-Stratonovitch transformation. Unfortunately 
such a transformation breaks planarity of the diagrams and is not suitable for our purposes. 
However, another Hubbard-Stratonovitch transformation exists for the particular values 
n = ±2 (besides n = 1, of course) which preserves planarity, as will be explained later. 

We can also define correlation functions in the model. There is only one 2-point 
function, 
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where a is fixed, 1 < a < n. Here the brackets (•) refer to the normahzed average with the 
exponentiai weight of (2.1). Due to 0(n)-invariance, there are two independent 4-point 
functions. We consider connected 4-point functions only; we choose 
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(2.5a) 
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where a and b are fixed and distinct. They have the following interpretation: ^(g) is 
the generating function of the number of colored alternating tangle diagrams of type i; 
a diagram with four external legs is of type 1 if the external strings come in and out in 
diagonally opposite corners; it is of type 2 if the external strings come in and out in the 
same upper/lower half-plane; see Fig. 2. Note that any diagram with four external legs is 
either of type 1, or of type 2, or the image of a type 2 diagram under a tt/2 rotation. 









Fig. 2: The two types of tangle diagrams. 



We must now proceed as in [2]. Our goal is 1) to count only diagrams which are prime 
and reduced, and 2) to count as a single contribution diagrams which correspond to the 
same link. First it is necessary to take care of non-prime and non-reduced diagrams. This 
is achieved by introducing an extra parameter t in the action: 



zW(n,t,g)= l[dM a e 



Ntv(-I ELi M a + f Ea tb =i M a M h M a M h ) 



(2.6) 



a=l 



and choosing t as a function of g in such a way that 

G(t(g),g) = l. (2.7) 

We have the obvious scaling property G(t,g) = jG(l, g/t 2 ), which means that given the 
two-point function of the original model G(l, g) = G(g), t(g) is the solution of the equation: 

t(g)=G(g/t 2 (g)) . (2.8) 

We have similar scaling properties for the higher correlation functions; in particular, 

r i (t,g) = ±r i (i, g /t 2 ). 




Fig. 3: The flype of a tangle 

The next step is to remove the over-counting of links due to the fact that several 
diagrams can correspond to a single link. According to the Tait flyping conjecture (proven 
in [8] ) , two reduced alternating diagrams are equivalent if and only if they are related by a 
sequence of flypes (Fig. 3). In order to take into account the flyping equivalence, we must 
now consider, as advocated in [9], the most general (9(n)-invariant model with quartic 
interaction 

**>(»,«,»,») = ffi dM.."" H E - Ml + * ^ MM > 2 + * M « 

J a=l 

(2.9) 

There is a new type of vertex which allows loops to "avoid" each other. The appearance 
of two types of vertices can be understood as follows: the unwanted (i.e. overcounted) 
tangle diagrams due to the flyping equivalence can be of either type 1 or type 2, and 
we must introduce "counterterms" of both types to cancel them. To pursue the analogy 
with renormalization, we can rephrase this by saying that it is only the "renormalized" 
coupling constants which must be of the form (g,0), but in order to reach this point we 
must consider more general "bare" coupling constants (gi,g2)- 
We impose again the condition 



G(t(g 1 ,g 2 ),gi,g 2 ) = 1 



(2.10) 



From now on t will be assumed to be fixed by (2.10), and we shall be concerned with 
the correlation functions = (f7i , #2), gi, #2), which are the generating functions of 
diagrams of type i with four external legs and two types of vertices weighted by gi and #2 
respectively 

As in [2] , we make use of the concepts of two-particle irreducibility. Following the lan- 
guage of field theory, we recall that a four-legged diagram is two-particle-irreducible (2PI) 
if cutting any two distinct propagators leaves it connected; otherwise it is two-particle- 
reducible (2PR). Also, we shall consider skeleton diagrams, which must be "dressed" to 
recover ordinary diagrams. This will be implicit in the following and we refer the reader 
to [2] for details. 

Let us define Di (resp. Hi, Vi) to be the generating functions of 2PI (resp. 2PI in 
the horizontal, vertical channel) four-legged diagrams of type i. Note that Hi = V±, but 
H2 7^ V2, since the defining property of diagrams of type 2 is not invariant by rotation of 
tv/2. 

By decomposing a general diagram I\ according to the number of times it is reducible 
in the horizontal channel (Fig. 4), we find the following formulae: 




(2.11a) 
(2.11b) 



r 2 H 2 H l H l H 2 H 2 

Fig. 4: Decomposition of a tangle diagram in the horizontal channel. 
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This can be simplified by introducing the combinations At = A ± A and similarly 
for # ± . We find: 

r ± = . (2.110 

Note that these equations are independent of n. 

It is not as obvious how to decompose diagrams of type 2 using the vertical channel. 
It is simpler to proceed as follows: define A to be generating function of diagrams with 
four external legs such that the color of the two left outgoing strings is free, whereas the 
color of the two right outgoing strings is fixed and equal; and similarly H 0: Vb, -Do- We 
have the formulae 

r = (n + l)r 2 + A (2.12a) 
H = H 2 + nV 2 + H x . (2.12b) 

We can now proceed to decompose A in the horizontal channel. It is easy to convince 
oneself that the simple formula 

ro = j\- (2.13) 
1 — Ho 

holds. 

The three equations (2.11) and (2.13) determine Hi = Vi, H 2 and V 2 as functions of 
A and F 2 } More explicitly, one should invert them to 

Ho,± = (2.14) 
1 + J- o,± 

and take appropriate linear combinations. Furthermore, noting that any diagram with 
four external legs and four-legged vertices must be irreducible in one of the two channels, 
we have 

D i = H i + V i -r i , (2.15) 

which means that we have managed to express the A in terms of the A- 

The flype equivalence does not affect 2PI skeleton diagrams; this means in practive 
that the expression of D\ = Di — gi as a function of the A is left unchanged by the 
removal of the overcounting of flype equivalent diagrams. Note that the Dj[A> A] have a 
simple expression given by Eqs. (2.14)-(2.15), but the <7i[A> A] are non-trivial functions 



1 If one sets n = 0, two equations become identical and one has to include the derivative of 
(2.13) with respect to n: A = V 2 /(l - (Hi + H 2 )f. 
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which can only be obtained by actually solving the matrix model (2.9) (or some equivalent 
procedure). For example, perturbatively we find: (Fig. 5) 

d[ = n rf + 8rfr 2 + {An + 4)rlr% + 24/frl + isrfr 2 + o( g 9 ) (2.16a) 
d' 2 = nrtr 2 + 2rl + i6r-frf + (20 + 8n)/?r| + (14 + 6n)rfr 2 + 3rf + o( g 9 ) 

(2.16b) 

where the order of truncation is dictated by the rule: 2 i~i ~ g, F 2 ~ g 2 . 




Fig. 5: First terms in the perturbative expansion of D[ and D' 2 . 

The data -D^-Ti,/^] is all we need from the generalized matrix model (2.9), and we 
can now come back to the problem of tangle diagrams, with a single coupling constant g 
weighting a simple crossing. 

We want to redo the decompositions of a general diagram in the horizontal channel, 
but this time taking into account the flyping equivalence. This forces us to distinguish 
between simple crossings and other diagrams, that is to use the primed objects defined by 
r[ = i\ — g, H[ = Hi — g, V{ = V\ — g. Here, as in [2], the tilde on one of the symbols T, 
H or V denotes generating functions of flype equivalence classes of 2PI skeleton diagrams. 

2 This rule can be justified as follows: at leading order, the bare coupling constants can be 
replaced with their renormalized values: g\ = g, gi = 0. One then finds 7i = g + 0(g 2 ), 
r 2 = g 2 + 0{g 3 ). 
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We find: 




(2.17a) 



r! 2 = g r 1 + \[ — ^ ' "\ + —„ "\ I . (2.17b) 

2 I 1-(H 2 + H[) 1-(H 2 -H[) 1 



Introducing H' ± = H' 2 ± Hi , we rewrite this 

(l T 2)f ± = ±2+— I* . (2.18) 
1 — H, 



Similarly, if H' Q = H — g, we find 



(l- g )f = g+-^k- . (2.19) 



We invert these relations to 



± l + (l T g)f T T9 

fi, = (1 - g)f Q - g 
l + (l- g )r -g 



and express the 2PI functions 



(2.20a) 
(2.20b) 



D[ = H[ + V{ - f[ = H' + - H'_ - P 1 - g (2.21a) 

D' 2 = H 2 + % -f^ = ±(H' + + H'_) + \H' - H' + ) - f 2 . (2.21b) 

Finally, the generating functions r^g) are given by the implicit equations 

b' i {g) = D' i [r 1 {g),r 2 {g)] (2.22) 

where D' i [ri,r 2 ] comes from the solution of the matrix model, and D'^g) is given by 
Eqs. (2.20)-(2.21) with f< = f^g). 

We can solve the implicit equation (2.22) perturbatively using (2.16); the result is 
summarized in Tab. 1. 





g 9 2 g 3 


9 A 


9 5 


9 6 


9 7 


9 8 


A 


1 2 


2 


6 + 3n 


30 + 2n 


62 + 40n + 2n 2 


382 + 106n + 2n 2 


f 2 


1 1 


3 + 71 


9 + n 


21 + lln + n 2 


101 + 32n + ?i 2 


346 + 153n + 24n 2 + n 3 



Tab. 1: Table of the number of prime alternating tangles up to 8 crossings. 
The power of n indicates the number of closed loops, i.e. the number of 
connected components besides the two outgoing strings. 
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These results are compatible with the exact solution at n = 1 [3], and they will give 
us a non-trivial check of our new solution at n = 2. 



3. The n = 2 case: a model of oriented links 

Let us now carry out explicitly the procedure outlined in the previous section, for 
a value of n corresponding to a solved matrix model. The case n = 1 has already been 
considered in [2]; let us now set n = 2. The partition function 

Z = J dMidM 2 e^ tr H( M i + M ') + ^(Mf + Mi) + f(M 1 M 2 ) 2 + g 2 M?Mi) 

(3.1) 

is conveniently rewritten in terms of a complex matrix X = * / 1 (Mi + iM<i): 



Z = J dXdXt e Nt * (~ XXt + b ° X2X ' 2 + ¥o{XX^) (3 2) 

where we have absorbed t in the coupling constants: bo = b/t 2 , Co = c/t 2 , with b = g\ + g% 
and c = 2<72- It is clear that the partition function has been left unchanged by the 
transformation, but the individual diagrams are different. The new Feynman rules are 
depicted on Fig. 6 a). If we simply set b = g, c = this model describes oriented links, 
each closed loop having now two possible orientations instead of two colors. 

A e i9 ' 2 -------f > 



V 

a) b) 

+ie/2 



/\ 

Fig. 6: Feynman rules corresponding to a) (3.2) and b) (3.3). 

We recognize in (3.2) the partition function of the six-vertex model on random dy- 
namical lattices, which has recently been solved [4,5]. In the 6 vertex formulation it is 
particularly clear that there is a Hubbard-Stratonovitch transformation which preserves 
planarity; 3 introducing the notation c = 26 cos 9, with 9 G [0, 7r[ in the regime of interest 
to us, we have 



f , Ntr (-XX* - ^-A 2 + A(XX^e l6 / 2 + X^Xe 

/ dALYdX f e V 2b ° v 



■id/2- 



(3.3) 



The same remark applies to the corresponding fermionic model, i.e. n = —2. 
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where A is a hermitean matrix. (Here and in the following integral (3.4), we omit an 
overall constant factor, irrelevant in the computation of correlation functions.) We have 
rewritten our model as a model of oriented loops which are not the intersecting loops we 
started from, since they can only avoid each other, see Fig. 6 b). This may seem unnatural, 
but is what allows to solve exactly the model, since we can now integrate over the original 
matrix X, shift the matrix A: 

Z = J dAdet - 1 (>/ 2 <g) A + A <g) e"^/ 2 ) e ~ N ?k tr ( A ~ 2cos(*/2)) (3.4) 

and do a saddle point analysis of the eigenvalues of A. This is the basis of the analytic 
solutions [4,5] of the model. 

It is natural to redefine the two independent 4-point functions to be 

n = ^tr(X 2 X^ (3.5a) 
r c =^tr(XX^ (3.5b) 

i.e. they are characterized by the position of the ingoing/outgoing arrows on the external 
legs. They are related to the 4-point functions defined earlier by: r b = Ti + T 2: T c = 2r 2 , 
as is clear diagrammatically. 

Let us summarize the combinatorial relations of the previous section in the case n = 2. 
First, the relations which determine the corresponding 2PI correlation functions D b and 
D c can be obtained either by direct diagrammatic arguments or by setting n = 2 in the 
general formulae (2.11)-(2.13). We find: 

A = r ^ r (3.6a) 

D' b = H b + V b -r b -b (3.6c) 
D' c = H c + V c - r c - c . (3.6d) 



Similarly, the modified relations which take into account the flyping equivalence read: (we 
have reintroduced the renormalized coupling constant g corresponding to a simple crossing, 
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and H' b = H b — g, etc) 

f b = g(l + f b ) + -?L- (3.7a) 
1 ~ H b 

(1 =F g) (f c ± P b ) = ±g + Vc± V K (3.7b) 

D' b = H b + V b -r b -g (3.7c) 
D' c = H c + V c -r c . (3.7d) 

We shall now use the solution of the matrix model (3.2) [4,5] to show how to extract 
the functions t perturbatively at an arbitrary order, and to find their large order behavior. 

3.1. Perturbative expansion of the off-critical solution. 

In [5] , the large N saddle point density of eigenvalues of A is explicitly constructed in 
terms of elliptic functions. More precisely if we define the resolvent 

W(a) = lim ( — tr^— > (3.8) 

where the average is with respect to the measure in (3.4), and 

J(a) = i \w(iae- ie ' 2 ) - W(-iae i6 / 2 )} + 3-— a 2 - * - a (3.9) 

v ; L v ; v >\ 2sin^6 4cos 2 (^/2)6 

then we have the following expression for J using an elliptic parametrization u: 

J = A + B^—^ (3.10b) 

sn 2 (-u--u oc ) 

a = a — (3.10a) 

Hiuoo - u) 

where H is the Jacobi theta function and A, B, ao are constants which depend on 6 
and 9. 

If we consider the small 6, c perturbative expansion of the correlation functions of 
the model, we are in the region where the elliptic nome q is close to zero, and the elliptic 
functions can be expressed to a given order in q in terms of trigonometric functions. We 
can therefore write every quantity as a power series in q with coefficients dependent on 9. 
From the 1/a term in the a — > oo expansion of J(a) we can extract W\ = liniAr^oo (-^trA), 
and from there the two-point function using the formula 

G =^G^r w '0' (3 - n) 

li 



We can then go back to the free energy by integrating once (at fixed 9) 

-b (7—1 

d6o ^T (3 ' 12) 

and differentiate again to get the 4-point functions: 

d 

H = ~QQ^\b fixed • (3.13) 

The rescaling which allows to set G = 1 simply amounts to writing that the coupling 
constant b is 

b = b G 2 . (3.14) 

Then we have 

rk= (G-w + Hm te_ 1 (3i5a) 
r ' = -FiL- 2 ' (315b > 

We must apply formulae (3.6) to evaluate D' b and D' c . Finally, we must consider q and 
as power series in the renormalized coupling constant g and solve perturbatively in g the 
equations (3.7). 

This can be programmed on a computer using for example Mathematica™. At the 
first non-trivial orders we find: 

b = q 2 - 6(1 + 2 cos % 4 + 0(q 6 ) 

G = 1 + 2(1 + 2 cos 6)q 2 + 0(q 4 ) 

r b = q 2 -q 4 + 0(q 6 ) 

r c = 2 cos 9 q 2 + 2(1 - cos % 4 + 0(q 6 ) 

D' b = (6 + 12 cos9 + 4 cos(29) + 4 cos(3e))q 10 + 0(q 12 ) 

D' c = (8 + 24cos^ + 8cos(2^) + 10cos(3#) + 2cos(5e))q 10 + 0(q 12 ) . 
The next step is to solve perturbatively the Eqs. (3.7). We can theoretically proceed to 
an arbitrary order in q and therefore in g. In Tab. 2 we show the results obtained after 
eight hours on a Sun work-station. As a check, note that g\ = (1 — cos 0)b = g + 0(g 4 ) 
and <72 = cos 6 b = —g 3 + 0(g 4 ) which is consistent with the fact that at order 3 there is 
exactly one diagram of type 2 (up to a it/2 rotation) which is overcounted (see Fig. 7). 
Also, these data are compatible with those of Tab. 1. 





Fig. 7 : First flype-equivalent tangle diagrams. 
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1 


1 


4 


13 


319 
6 
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2 


1941 
2 


13424 
3 


ooozu J 
40 


o44o 1 1 
8 


12 


ni orni nr 
zlooa lUo 

8 


112 


louyiJoiuyj 
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b 


1 




— 1 


—3 


—9 


—27 


— 103 


—411 


— 1838 


—8484 


—41000 


—202822 


—1027954 


—5300828 


—27768068 


—147347688 


c 






-2 


-2 


-6 


-18 


-74 


-314 


-1420 


-6696 


-32592 


-162628 


-828244 


-4290608 


-22553488 


-120042224 


9i 


1 






-2 


-6 


-18 


-66 


-254 


-1128 


-5136 


-24704 


-121508 


-613832 


-3155524 


-16491324 


-87326576 


92 






-1 


-1 


-3 


-9 


-37 


-157 


-710 


-3348 


-16296 


-81314 


-414122 


-2145304 


-11276744 


-60021112 


h 


1 


1 


3 


7 


23 


81 


319 


1358 


6132 


28916 


140852 


704020 


3592394 


18648750 


98217650 


523734140 






2 


2 


10 


22 


94 


338 


1512 


6700 


31944 


155200 


778168 


3972088 


20646268 


108809588 


580654216 


A 


1 




2 


2 


12 


34 


150 


602 


2782 


12944 


63252 


314936 


1606350 


8325616 


43812856 


233407032 


h 




1 


1 


5 


11 


47 


169 


756 


3350 


15972 


77600 


389084 


1986044 


10323134 


54404794 


290327108 



Tab. 2: Table of the number of oriented prime alternating tangles up to 16 
crossings, and related quantities. A# = 9 — tt/2. 



3.2. Solution at criticality and large order behavior of the correlation functions. 

In [4], the model (3.2) is analyzed in the critical regime, that is when one encounters the 
first singularity of the free energy (or of the correlation functions) by analytic continuation 
from the gaussian model. This singularity - or more precisely, the line of singularities in 
the (&o, Co) plane - is what determines the large order behavior of the correlation functions, 
and ultimately of the functions r. 

We shall now summarize some of the relevant results from [4]. If the parameter 9 is 
fixed and the coupling constant bo is increased, for any value of 9 G [0, 7r[ there will be a 
critical value &o(6>) for which the free energy becomes singular: 

1 tan(fl/2) 1 mfil 
bo{9) ~ 32 9/2 cos^/2) • (3 - 16) 

This forms a line of critical points, describing what is known in physics as the line of 
c = 1 conformal field theories - i.e. a free boson compactified on circles of varying radius 
- coupled to gravity, see for example [10] for a review and references. In [4], all correlation 
functions of the form lim/v^oo (jj^ rj ^ s ) are calculated on the critical line. It is easy to see 
that the explicit expressions of W\ = liniTv^oo (]^trA) and W2 = limjv^oo (-^trA 2 ) allow 
to extract the 2-point function and the free energy (still on the critical line); one finds 

By definition F*(9) = F(6q(#)), 9), which means that 

A 1 Ah* F) 

This gives us access to H = -§gF\ bo fixed at criticality. One finds: 

H *(9) = -29 +l(n 2 - 9 2 ) tan(0/2) - \n 2 \ + hn 2 - 9 2 ) cot(#/2) . (3.19) 
2 6 9 b 

One must then use again the formulae (3.14)-(3.15) and (3.6) to find the values of -T c , 
D' b , D' c on the critical line as a function of 9, and finally solve (3.7) for g and 9. This is a 
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set of 2 complicated coupled equations, 4 and we can only solve them numerically. We find 
the following numerical values: 

9 C = 1.60780446 . . . (3.20a) 
l/g c = 6.28329764... (3.20b) 

The value (3.20b) is non-universal, and only this explicit calculation could give us access 
to it. It is related to the leading behavior of the coefficients of the power series r b (g) 
and r c (g), i.e. of the numbers of oriented prime alternating tangles. On the other hand, 
the subleading behavior is universal; it is characteristic of c = 1 conformal field theories 
coupled to gravity, which exhibit zero "string susceptibility" with logarithmic corrections 
[4,5]. Therefore we can state that if /&(#) = J2lp9 P i then 

7 /~°°const g-Pp-^logp)- 1 (3.21) 

and similarly for P c . The constant l/g c = 6.28329764. . . is slightly less than the corre- 
sponding constant 16/(7r(7r — 4) 2 ) = 6.91167. . . [9] for oriented alternating link or tangle 
diagrams (without taking into account the flype equivalence); and should also be compared 
to the numbers obtained for n = 1: (101 + v / 2100T)/40 = 6.147930 ... and 27/4 = 6.75 
with and without the flype equivalence, respectively. The fact that the n = 1 and n = 2 
results are fairly close shows in particular that the entropy generated by tangles with large 
numbers of connected components is small. 

As a final note, we can make some slightly conjectural statements on the number of 
oriented prime alternating links. It is clear that if we close a tangle by pasting a simple 
crossing to its four external legs, we obtain a link; if the tangle had k closed loops, then 
the number of connected components of the resulting link is simply k + 2 for a tangle of 
type 1 and k + 1 for a tangle of type 2. Inversely, from a link with p crossings one can 
produce 4p tangles by removing one of its vertices and fixing the circular permutation of 
the four legs. Even though we cannot use this fact to correctly count links (one cannot 
make it a one-to-one correspondence), by assuming that most links have low symmetry, 
we have the approximate relation f p 7 P _i/p satisfied by the number f p of links with p 
crossings. Therefore we conjecture that 

//~°°const g-Pp-^logp)- 1 . (3.22) 

4 In fact, one of the two equations is quadratic in g, so that remains only one (transcendental) 
equation in 9. 
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